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We consider a reaction-diffusion model incorporating the reactions A — > 0, A — > 2A and 2A — > 3A. 
Depending on the relative rates for sexual and asexual reproduction of the quantity A, the model 
exhibits either a continuous or first-order absorbing phase transition to an extinct state. A tricritical 
point separates the two phase lines. As well as briefly examining this critical behavior in 2+1 
dimensions, we pay particular attention to the cluster geometry. We observe the different cluster 
structures that form at criticality for the three different types of critical behavior and show that 
there exists a linear relationship for the probability of survival against initial cluster size at the 
tricritical point only. 
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Being able to identify the important geometrical struc- 
tures is known to greatly facilitate the understanding of 
the underlying physical mechanism [1]. In this letter we 
describe how the tricritical point in a certain class of 
nonequilibrium population models is characterized by the 
unique structure of the spatial cluster. Nonequilibrium 
phase transitions are of great relevance to disciplines as 
wide-ranging as traffic flow [2] , chemical reactions and 
atmospheric studies [H . Of particular interest is the idea 
of universality, where different models display identical 
scaling functions and critical exponents close to the crit- 
ical point. By far the most well-known and studied uni- 
versality class out of equilibrium is directed percolation 
(DP). The robustness of DP led Janssen and Grassberger 
@, E| to the conjecture that all models with a scalar or- 
der parameter that exhibit a continuous phase transition 
from an active to a single absorbing state belong to the 
class. A universality class closely associated with DP is 
tricritical directed percolation (TDP), which has recently 
been studied both from steady-state [3] and dynamical Q 
simulations. 

The process of TDP incorporates higher-order terms 
than DP, with the multicritical behavior occurring if the 
lower-order reactions vanish on a coarse-grained scale Q . 
Liibeck [7[ examined tricritical behavior by adding the 
pair reaction 2A — > 3A to the contact process [10] while 
Grassberger [8fl used a generalisation of the Domany- 
Kinzel model [ll| in 2+1 dimensions. 

Recently, we examined a simple reaction-diffusion 
model with the reactions A — > <j) and 2 A — > 3 A [l2j], 
which exhibited a phase transition that is continuous, 
and DP in 1+1 dimensions and first first-order in higher 
dimensions. Although, from the mean field, the first- 
order phase transition was expected in all dimensions, 
the larger fluctuations in the (l+l)-dimensional case are 
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FIG. 1: Typical plots for sexual reproduction only (left) and 
different rates of both sexual and asexual reproduction (mid- 
dle and right). 



likely to destabilize the ordered phase [13], resulting in 
the observed DP transition. 

Following on from Liibeck [7], we introduce a lower- 
order proliferation reaction to our original model. By 
allowing both asexual and sexual reproduction, we have 
a phase space exhibiting first-order, continuous and 
tricritical behavior. Whilst the dynamics of such a 
system have been previously studied [7|, [8j , the geomet- 
rical structure of the population at the different types 
of critical behavior has so far been ignored. Sample 
snapshots of the population for sexual reproduction only 
and both asexual and sexual reproduction are shown in 
Fig. [TJ In this letter, we study this geometrical structure 
and show the remarkable result that the relationship be- 
tween initial cluster size and the survival probability up 
to some finite time t is linear at the tricritical point only. 

The Model: We have a d-dimensional square lat- 
tice of linear length L where each site is either occupied 
by a single particle or is empty. A site is chosen at 
random. The particle on an occupied site dies with 
probability pa, leaving the site empty. If the particle 
does not die, a nearest neighbor site is randomly chosen. 
If the neighboring site is empty the particle moves there 
and produces a new individual at the site that it has 
just left with probability k. If the chosen site is however 
occupied, the particle reproduces with probability 
producing a new particle on another randomly selected 



2 



neighboring site, conditional on that site being empty. 
A time step is defined as the number of lattice sites 
N = L d and periodic boundary conditions are used. 
We have the following reactions for a particle A for 
proliferation and annihilation respectively, 



A + A- 



3 A, A- 



2A and A 



0. (1) 



Assuming the particles are spaced homogeneously, the 
mean field equation for the density of active sites p(t) is 
given by: 



dp(t) 
dt 



= Pb (l-p d )p(t) 2 (l-p(t)) 

+ k(l- PA )p(t)(l-p(t))-p dP (t). (2) 



The first term considers the sexual reproduction term, 
the second asexual reproduction, and the final term death 
of an individual. Eq. [2] has three stationary states: 



po 0, 



(3) 
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For k > pb, P+ — » continuously as pd — > fc/(l + fc), 
indicative of a continuous phase transition with critical 
point 



Pd c 



1 + k 



(4) 



For k < however, we have a jump in p + from (pb 
k)/2p\y to zero, this time at the critical point 



Pd c 



(k+p h ) 2 



4p b + (A:+Pb) 2 ' 
Further, for k < pb : we have a region 

k . (k+p h ) 2 



< Pd < 



4p b + (&+Pb) 2 



(5) 



(6) 



where the survival of the population is dependent on the 
population density. In fact, we have extinction for 



p(t) < p-(k <p h ,Pd) 



(7) 



since then dp/dt < 0. For k < Pb we therefore have 
a first-order phase transition. The two lines of phase 
transitions meet at the point k = pb, defining the position 
of the tricritical point k*. At the mean field level then, we 
have a phase diagram as shown in Fig. [2j We note that in 
our earlier paper [12[ , we examined the case k = and so 
were restricted, at the mean field level, to the first-order 
regime only. 

From simulations we find that for 1+1 dimensions, con- 
tinuous phase transitions are observed for all k. In order 
to examine tricritical behavior we therefore proceed in 
2+1 dimensions. 
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FIG. 2: Phase diagram showing the continuous phase transi- 
tion occurring for k > ph and the first-order phase transition 
for k < ph. Eqns. JU|5} have been re-arranged in the figure 
to make k the dependent variable. 



To examine the critical behavior of our model, we 
study the dynamical behavior of the system after starting 
from a single- seed at the centre of a large lattice. The 
lattice is sufficiently large so that no individuals reach 
the boundary during the running of the simulation. For 
fixed pb = 0.5 we find the critical value p^ c for a given 
k. For the continuous phase transitions, at p^ = pd c , we 
expect power-law behavior 



n(t) 6c P 
P{t) 6c t~ ( 



(8) 
(9) 



for the population size and probability of survival re- 
spectively. For k sufficiently greater than fc*, we find 
DP values r] = 0.231 and 5' = 0.451 [14] as shown in 
the inset of Fig. [3j At k = k* we again expect power- 
law behavior because of the continuous phase transition, 
but with different values for the exponents 0, [l5[ . The 
best power-law is found for k = k* = 0.12(0) with 
Pd = Pd c = 0.150219(5) (see Fig. [3j), where the num- 
bers in the perenthesis show the uncertainty in the last 
figure. For more details of methods obtaining the tricrit- 
ical point, see Refs. [7|, l8|]. 

To find the values of pa c in the first-order regime 
we adopt a different approach, inspired by Lee and 
Kosterlitz [16[, due to the lack of power law behavior. 
For k < k* we examine histograms of population density 
against frequency for different values of pd- Due to the 
phase-coexistence at the critical point, we observe a 
double-peaked structure in the histogram for p^ close to 
Pd c , where the peaks are at equal heights for pd = Pd c • 

Cluster Geometry: The difference between the 

mean field and simulation results is due to both the 
neglect of the fluctuations in Eqn. ([2]) and the false 
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FIG. 3: Plots showing the power- law behavior for n(t) and 
P(t) (bottom) outlined in Eqns. ([HE]). The main plot is 
with k = k* = 0.12 and p d = p dc = 0.1502195. We find 
rj — —0.487(5) and 5' — 1.36(5). The inset shows (from top 
to bottom) k = 0.5, 0.75 and 1. The values of pd c are found 
to be 0.30302(0), 0.36765(0) and 0.41626(0) respectively. The 
hashed lines show the power-law with DP values. 

a) 



b) 



FIG. 4: Typical plots of the population structure at the criti- 
cal point for a) (from left to right) k = 0, 0.05, 0.5 and 1 and 
b) k = k* = 0.12. The simulations began from an initially 
fully occupied 128 x 128 lattice and ran until the population 
size reached 400 for the first time. 

assumption of a homogenous distribution. In this letter, 
we are particularly interested in the heterogeneous 
distribution and in particular the geometrical structure 
of the clusters at criticality and how this changes with 
k. From Fig. [Ha), we see that for fccF, after starting 
from a fully occupied lattice, the population becomes 
densely clustered. This is certainly not unexpected 
since, for A: = 0, reproduction is only possible if an 
individual finds another. In contrast, in the DP region, 
the population is much more dispersed. For k ^> fc*, an 
individual has a greater chance of reproduction if it is 
surrounded by empty sites and so being heavily clustered 
is a disadvantage. For k ^ k* , the plots shown in Fig. 
2] a) are very typical plots, whereas at the trier it ical 
point, the plots are very varied. For k = k* , Fig. [H b) 
shows that the population has little preference as to its 
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FIG. 5: Plot of survival probability (for t ma x = 200), against 
initial cluster size for (from bottom to top) k = 0, 0.12 and 1. 
Relatively small cluster sizes were used since clearly P s — > 1 
as cluster size — > oo. 



structure. There are often areas where the population is 
heavily clustered as well as those where the structure is 
DP-like. 

In considering the geometrical make-up of the clusters, 
it is interesting to examine how the survival probability 
of a cluster depends on its size for different values of k. To 
examine this, we run simulations from an initially fully- 
occupied lattice up to the point when a certain popula- 
tion is reached. We ensure that this population is suffi- 
ciently small so that the resulting density does not 'force' 
the population to take on certain cluster sizes. Once this 
population has been reached, we make a copy of each 
cluster and place each one in the centre of a sufficiently 
large lattice. The survival probability P s up to some time 
t max of the population is then measured for each individ- 
ual initial cluster. The results for each cluster size are 
shown in Fig. [5l 

We see that the gradient of P s with respect to cluster 
size increases for k = whilst it decreases for k = 1. 
Consistent with the plots from Fig. 0J we have that ad- 
ditional individuals give an increasing advantage to the 
survival probability for k = but a negative effect for 
k = 1. Due to the curvature of the plots for k = 1 and 
k = 0, we would expect a linear relationship for some 
< k < 1. Surprisingly this occurs at the tricritical 
point k = k* = 0.12. In order to understand why this 
linear behavior might occur at the tricritical point, we 
plot in Fig. [6] the average number of births and deaths 
that occur per individual per time step, B\(n) and V\(n) 
respectively, at the critical point for different k. As ex- 
pected, for all values of /c, T>i(n) remains constant for all 
population sizes. In contrast, an increasing population 
seems to have a positive influence on B\(n) for k = 
whilst a negative influence for k = 1. At k = for suf- 
ficiently large population sizes, B[(n) is independent of 
the population size. In fact, for n > 40, B\{n) w V\(n). 
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FIG. 6: The number of births (solid line) and deaths (hashed 
line) per individual per time step for (from top to bottom) 
k = 1, 0.12 and 0. 



We examine now whether it is this equality in the num- 
ber of births and deaths that results in the observed lin- 
ear relationship between initial cluster size and P s . We 
simplify our model by considering the macroscopic pop- 
ulation size n(t) only and ignoring the microscopic in- 
dividuals. Analogous to a random walk, the population 
increases by one with probability £>b> decreases by one 
with probability pa and stays the same with probability 
1 — pb — Pd- We again examine the probability that the 
population survives up to some time t max where a time 
step at time t is defined as n(t) updates. We note that 
this is the same definition of a time step that we use when 
running simulations that begin with a single-seed. In the 
simulations, it corresponds to approximately one update 



per individual. 

With p\^ = pa, we find that the probability of sur- 
vival increases linearly with initial population size. This 
finding seems to indicate the equality in the number of 
births/deaths per individual per time step results in the 
observed linear behavior at the tricritical point. 

Conclusion: We have examined the critical behavior 
in our model which exhibited continuous, first-order and 
tricritical phase transitions. We noted the key geomet- 
rical differences in the cluster structure for each type of 
transition caused by the different values of the rate of 
asexual reproduction k. Of particular importance was 
the surprising result of the linear increase in probability 
of survival with cluster size. We hypothesised that this 
was due to the number of births and deaths per time step 
being equal at the tricritical point. 

The lack of sensitivity of this linear behavior to the 
parameters do not make this an effective method to ob- 
tain the position of the tricritical point. To do this, we 
believe seeking the power law behavior for n(t) and P(t) 
to be the best approach. 

One area of future study is the discrepency between 
our values for n and 8' at the tricritical point with Grass- 
berger's who found rj = -0.353(9) and 8' = 1.218(7) @. 
Indeed, he had disagreement between his values for /?, 
v± and <j) and Liibeck's [7]. Universality at the tricritical 
point seems to have been assumed but the differences in 
these critical exponents appear to dispute that. Further 
study needs to be done to clarify the situation. 

We are indebted to Sven Liibeck whose correspondence 
to us about our first paper on this subject led us to 
this research. All computer simulations were carried out 
on the Imperial College London's HPC. Alastair Win- 
dus would also like to thank EPSRC for his Ph.D. stu- 
dentship. 
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